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Fig.  1.  Spin  echo  images  produced  at  successive TR  times  with  shortTE. 


The  images  were  taken  from  slices  of  4  mm  thickness 
and  their  sizes  are  256  x  256. 

The  pixel  values  in  16  images  were  plotted  against  TR 
values  as  shown  In  Fig.  2. 


Fig.  3.  A  sample  dataset  of  estimated  parameters  and  the  best  fit. 


The  fitting  algorithm  was  applied  to  all  pixels  and  (Ti, 
PD)  values  were  recorded  for  each  coordinate.  By  using 
these  estimated  values  the  calculated  Ti  and  calculated  PD 
image  (see  Fig.  4)  are  created. 


Fig.  4.  Calculated  Ti  and  calculated  PD  images. 

The  image  on  the  left  Is  the  trueTi  Image,  next  to  it  is  the 
true  PD  image. 


B.  Segmentation  of  Weighted  and  Calculated  Images 


Fig.  2.  Pixel  value  versus  TR. 

At  the  next  step,  Levenberg-M  arquardt  M  ethod  is  applied 
to  the  data  and  T i  and  PD  values  were  estimated. 

A  sample  data  together  with  the  fitted  exponential  curve 
and  the  estimated  parameters,  PI  and  P2  (they  stand  for  PD 
andTi)  is  shown  in  Fig.  3. 


In  my  research  1  used  maximum  likelihood  method  to 
choose  the  class  that  each  pixel  of  our  Image  belongs.  First, 
two  scatter  diagrams  are  produced  where  each  data  point  is  a 
function  of  Ti  and  PD.  One  of  these  scatter  diagrams  is 
achieved  by  using  the  Ti-weighted  (TR1=400  ms.)  and  PD- 
weighted  (TR16=1150  ms.)  images  (Fig.  5). 


Fig.  5.  Scatter  diagram  of  pixels  chosen  from  Tr weighted  and  PD-weighted 

images. 


Fig.  7.  The  multivariate  Gaussian  density  function  of  the  weighted  dataset. 


Blue  dataset  is  CSF,  red  dataset  is  GM ,  and  green  dataset 
is  WM.  The  training  dataset  that  wiil  be  used  in  weighted 
image  segmentation  is  created  by  choosing  some  coordinates 
on  the  brain  si  ices.  The  same  coordinates  are  aiso  necessary 
to  set  up  the  training  dataset  that  wiil  be  used  in  caiculated 
image  segmentation.  That  training  dataset  is  shown  in  Fig.  6. 


Fig.  6.  Scatter  diagram  of  pixels  chosen  from  trueTi  and  true  PD  images. 

Next  the  multivariate  Gaussian  density  functions  (Fig.  7 
and  Fig.  8)  of  both  the  weighted  and  the  calculated  datasets 
are  found  by  using  (1). 


V(2;r)"det(C,J 


(1) 


where  x  and  are  each  2  by  1  vectors  with  components  X|< 
and  E{XJ,  respectively,  for  l<k<2  and  is  the 

associated  2  by  2  covariance  matrix.  The  symboi  det(C^) 
denotes  the  determinant  of  the  covariance  matrix. 
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Multivariate  Gaussian  Density  Curves  ofTissues  in  Calculated  Images 


Calculated  Proton  Density 


Fig.  8.  The  multivariate  Gaussian  density  function  of  the  calculated  dataset. 


The  segmentation  aigorithm  based  on  the  maximum 
likelihood  ciassification  method  that  utiiizes  the  muitivariate 
Gaussian  density  functions  is  appiied  to  the  4  human  head 
slices,  one  of  which  is  iliustrated  in  Fig.  9  and  Fig.  10. 
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Fig.  9.  Segmented  slice  by  using  weighted  dataset. 


IV. Conclusion 


Fig.  10.  Segmented  slice  by  using  calculated  dataset. 


III.  Results 


To  evaluate  the  correctness  of  the  segmentation 
algorithm  the  confusion  matrix  method  is  used.  The  test 
pixels  are  chosen  from  the  coordinates  that  are  not  used  when 
setting  up  the  training  datasets.  To  be  certain  at  doing  this, 
the  training  datasets  are  constituted  from  the  tissues  on  the 
left  lobe  of  the  brain,  and  the  test  pixels  are  chosen  from  the 
tissues  on  the  right  lobe  of  the  brain.  The  segmentation  done 
with  the  calculated  datasets  gave  more  satisfactory  results 
than  the  segmentation  done  with  the  weighted  datasets. 

In  this  study,  development  of  a  segmentation  algorithm 
by  using  calculated  Ti  and  calculated  PD  images  is  aimed. 
For  this  purpose  the  true  T i  and  true  PD  images  are  produced 
by  finding  the  best  exponential  curves  fitting  to  the  images 
taken  from  an  M  R  imager,  and  the  segmentation  is  achieved 
by  using  the  maximum  likelihood  method.  We  conclude  that 
the  segmentation  done  with  the  calculated  datasets  gave  very 
satisfactory  results  and  that  it  will  be  helpful  in  detecting  the 
tumors  in  the  brain. 


To  evaluate  the  reliability  of  the  segmentation  algorithm, 
400  test  pixels  are  chosen  from  each  tissue  (CSF,  WM ,  and 
GM).  The  newly  chosen  pixels  (test  pixels)  are  different 
from  the  ones  that  constituted  the  training  dataset.  They  are 
used  as  input  to  the  segmentation  algorithm  and  the  results  of 
segmentation  are  arranged  in  confusion  matrices  as  shown  in 
Table  1  and  Table  2.  The  tissues  from  which  the  test  pixels 
are  chosen  are  placed  at  the  first  column  of  these  tables.  The 
first  row  shows  the  distribution  of  segmentation  among  three 
tissues.  For  example,  in  the  second  row,  348  test  pixels  are 
classified  as  CSF,  52  test  pixels  are  classified  as  GM ,  and  0 
test  pixel  is  classified  as  WM  from  400  test  pixels  chosen 
from  CSF. 


Table i 


SEGMENTATION  RESULTS  OF  400 TEST  PIXELS 
(BY  USING  WEIGHTED  DATASET). 


Result  of  segmentation 

Tissue 

CSF 

GM 

WM 

Tissues 

from 

which  the 

CSF 

348 

52 

0 

GM 

0 

356 

44 

pixels  are 
chosen 

WM 

0 

40 

360 
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Table  ii 


SEGMENTATION  RESULTS  OF  400 TEST  PIXELS 
(BY  USING  CALCULATED  DATASET). 


Result  of  segmentation 

Tissue 

CSF 

GM 

WM 

Tissues 

from 

which  the 

CSF 

360 

40 

0 

GM 

0 

388 

12 

pixels  are 
chosen 

WM 

0 

8 

392 

